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C^> Abstract 

The inclusive production of jets in the central region of rapidity is studied 
in fcT-factorization at next-to-leading order (NLO) in QCD perturba- 
' tion theory. Calculations are performed in the Regge limit making use 

Ph! of the NLO BFKL results. A jet cone definition is introduced and a 

proper phase-space separation into multi-Regge and quasi-multi-Regge 
^ ' kinematic regions is carried out. Two situations are discussed: scatter- 

ed , ing of highly virtual photons, which requires a symmetric energy scale to 

separate the impact factors from the gluon Green's function, and hadron- 
hadron collisions, where a non-symmetric scale choice is needed. 

X" 

1 Introduction 

The understanding of the physics behind jet production in the context of per- 
turbative QCD is an essential ingredient in phenomenological studies at present 
and future colliders. At high energies the theoretical study of multijet events 
becomes an increasingly important task. In the context of collinear factorization 
the calculation of multijet production is complicated because of the large num- 
ber of contributing diagrams. There is, however, a region of phase space where 
it is indeed possible to describe the production of a large number of jets: the 
Regge asymptotics (small-x region) of scattering amplitudes. This corresponds 
to the case where the center-of-mass energy in the process under study, s, can 
be considered asymptotically larger than any other participating scale. In this 
limit the dominating diagrams are those with gluons being exchanged in the 
i-channel. A perturbative analysis of these diagrams shows that it is possible 
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to resum contributions of the form (a s In s)™ to all orders, with a s being the 
coupling constant for the strong interaction. This can be achieved by means of 
the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation pQ. 

An essential ingredient in the BFKL approach is the concept of the Reggeized 
gluon or Reggeon. In Regge asymptotics colour octet exchange can be effectively 
described by a i-channel gluon with its propagator being modified by a mul- 
tiplicative factor depending on a power of s. This power, also known as gluon 
Regge trajectory, depends on the transverse momenta of the gluon and is not 
infrared finite. However, when real emissions are included using gauge invariant 
Reggeon-Reggeon-gluon couplings, the divergences cancel out. It is then possi- 
ble to describe scattering amplitudes with any number of particles (jets) in the 
final state. The (a s Ins)™ resummation is known as leading-order (LO) approx- 
imation and provides a simple picture of the underlying physics. Nevertheless 
it is not free of drawbacks, the main two being that, at LO, both a s and the 
factor scaling the energy s in the resummed logarithms, so, are free parameters 
not determined by the theory. These limitations can be removed if the accuracy 
in the calculation is increased, and next-to-leading (NLO) terms of the form 
a s (a s Ins)™ are taken into account [2J. When this is done, diagrams contribut- 
ing to the running of the coupling have to be included, and also So is not longer 
undetermined. As an example, in the context of Mueller-Navelet jets, the in- 
troduction of NLO effects in the kernel has been recently shown to have a large 
phenomenological impact, in particular, for azimuthal angle decorrelations |3]. 

At LO every Reggeon-Reggeon-gluon vertex corresponds to one single gluon 
emission, and the produced gluon can form a single jet. At NLO the situation 
is more complicated since the emission vertex also contains Reggeon-Reggeon- 
gluon-gluon and Reggeon-Reggeon-quark-antiquark contributions. In the present 
work we are interested in the description of the inclusive production of one jet 
in the BFKL formalism at NLO. This means that the relevant events are those 
with only one jet produced in the central rapidity region of the detector. In 
order to find the probability of production of a single jet it is necessary to intro- 
duce a jet definition in the emission vertex. This is simple at LO, but at NLO we 
should carefully study the possibility of a double emission in the same region of 
rapidity, leading to the production of one or two jets. This will be the main goal 
of the present paper. Our aim is to clearly separate the different contributions 
to the cross section, and to explain in detail which scales are relevant. Particu- 
lar attention is given to the separation of multi-Regge and quasi-multi-Regge 
kinematics. An earlier analysis has been presented in Ref. [4]. We have inde- 
pendently repeated these calculations, and we have found several discrepancies 
which will be explained in the text. 

Our analysis will be done in two different cases: inclusive jet production in 
the scattering of two photons with large and similar virtualities, and in hadron- 
hadron collisions. In the former case the cross section has a factorized form 
in terms of the photon impact factors and of the gluon Green's function which 
is valid in the Regge limit. In the latter case, since the momentum scale of 
the hadron is substantially lower than the typical &t entering the production 
vertex, the gluon Green's function for hadron-hadron collisions has a slightly 
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different BFKL kernel which, in particular, also incorporates some fc^-evolution 
from the nonperturbative, and model dependent, proton impact factor to the 
perturbative jet production vertex. We provide analytic formulae for these two 
processes, and the numerical analysis is left for a future publication. 

In the case of hadron-hadron scattering, our cross section formulae contains 
an unintegrated gluon density which, in addition to the usual dependence on 
the longitudinal momentum fraction typical of collinear factorization, carries an 
explicit dependence on the transverse momentum kx- This scheme is known as 
fcT _ factorization. So far, no systematic attempt has been made to generalize 
this framework beyond LO accuracy. In the small-a: region, where this type 
of factorization has attracted particular interest, the BFKL framework offers 
the possibility to formulate, in a systematic way, the generalization of the kx~ 
factorization to NLO. We therefore interpret our analysis also as a contribution 
to the more general question of how to formulate the unintegrated gluon density 
and the fey-factorization scheme at NLO: our results can be considered as the 
small-x limit of a more general formulation. 

After this short introduction, in Section 2 we define, closely following Ref. [5] , 
our notations for the description of a general cross section in the BFKL ap- 
proach. We also introduce multi-Regge kinematics (MRK) and the iterative 
structure of the cross sections at LO. In Section 3 we describe the basic ele- 
ments contributing at NLO. The linearity of the BFKL equation remains the 
same while the emission kernel now has several pieces such as virtual contribu- 
tions to one gluon emission and double emissions. We describe them in some 
detail, including a procedure to avoid double counting when the MRK is sep- 
arated from the quasi-multi-Regge kinematics (QMRK). The discussion of in- 
clusive jet production starts with a LO description in Section 4. Following this 
introductory part, we present, in Section 5, a definition of the NLO jet vertex. 
We separate the different regions of phase space in such a way that the can- 
cellation of infrared divergences is explicit for the two cases above-mentioned: 
inclusive jet production in 7*7* and in hadron-hadron interactions. We will 
also discuss the definition of a NLO unintegrated gluon distribution valid in the 
small-x regime. To close we study in Section 6 the role of the scale separating 
MRK from QMRK and show how, even with the jet definition, it is possible to 
prove that the dependence on this scale is power suppressed. Finally, we draw 
our Conclusions and suggest future lines of research. 

2 General structure of BFKL cross sections 

For the sake of clarity, in the present section we introduce the notation we 
will follow in the rest of this study. BFKL cross sections present a factorized 
structure in terms of a universal Green's function, which carries the dependence 
on s, and impact factors, which have to be calculated for each process of interest. 
This factorization remains unchanged in the transition from LO to NLO. We 
start by defining our normalizations at LO in the following. 

Lets consider the case of the total cross section ctab in the scattering of two 
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particles A and B. It is convenient to work with the Mellin transform 



f°° ds ( a\ ~ u 

J r (w,8 )= / — — CTAB, (1) 

Jso S \sqJ 

acting on the center-of-mass energy s. The dependence on the scaling fac- 
tor so belongs to the NLO approximation since the LO calculation is formally 
independent of sq. 

If we denote the matrix element for the transition A+B — > A+B+n produced 
particles with momenta ki (i = 1, ...,n) as ^g +n , and the corresponding 
element of phase space as #^g +n , then we can write 

1 °° f 

CTAB = Ys S / d$ AB+„I^AB+„| 2 - (2) 

As we mentioned in the Introduction we are interested in the Regge limit where 
s is asymptotically larger than any other scale in the scattering process. In this 
region the scattering amplitudes are dominated by the production of partons 
widely separated in rapidity from each other. This particular configuration of 
phase space is known as multi-Regge kinematics (MRK). In MRK produced par- 
ticles are strongly ordered in rapidity but there is no ordering of the transverse 
momenta which are only assumed not to be growing with energy. 

We fix our notation in Fig. [I] qi correspond to the momenta of those particles 
exchanged in the i-channel while the subenergies Sj-i,i = (ki-i+ki) 2 are related 
to the rapidity difference between consecutive s-channel partons. Euclidean 
two-dimensional transverse momenta are denoted in bold. For future discussion 
we use the Sudakov decomposition fcj = cxiPa 4- (3iPb + k i± for the momenta of 
emitted particles. 

In MRK the center of mass energy for the incoming external particles can 
be expressed in terms of the internal subenergies as 



n+1 

n 



n+1 



qx+i n -t==> ( 3 ) 



where we have used the fact that in Regge kinematics s is much larger than —t 
and, therefore, ao — (3 n +i — 1, M) ~ qf and k^ +1 ~ q n+1 . To write down the 
measure of phase space we use dimensional regularization with D = 4 + 2 e, i.e. 

i—l 5 \ / 

The matrix element A^+ n of Eq. ([2]) can be written in MRK in the factorized 
form 
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Figure 1: Notation for particle production in MRK. 

with Tp being the couplings of the Reggeon to the external particles, u>i = w(qf) 
the gluon Regge trajectory depending on the momentum carried by the Reggeon, 
and 7(^,^+1) the gauge invariant effective Reggeon-Reggeon-gluon vertices. 
At LO the scale sr is a free parameter. 

Gathering all these elements together it is possible to write the Mellin trans- 
form of Eq. (fTJ) as the sum 



71=0 



with the contributions from the emission of n s-channel gluons being 



(6) 



■r<"V,ao) 
(2tt) 2 - d 




n 2 



(7) 



The impact factors <£> p and the real emission kernel for Reggeon-Reggeon into 
a s-channcl gluon K, r can be written in terms of the square of the vertices T p 
and 7, respectively. The kernel /C r (qj,qj+i) is defined such that it includes 

one gluon propagator on each side: (q?q 2 + i) ■ The integration over Sj_ i,, in 
Eq. (0 takes place from a finite so to infinity. At LO terms of the form ui In kf 
or o^lnsfl can be neglected when the integrand is expanded in a s . Therefore, 
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at this accuracy, Eq. ([7]) gives 



jT(n)(, 



s J 



(27T) 



2-23 
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lD-2, 



2 0Ji 



■Mqi 



J^/C r (qi,q i+ i) 



i=l 



ga(qn+i) 

9 - 

q„+i 



(8) 



where the poles in the complex w-plane correspond to Reggeon propagators. 
This simple structure is a consequence of the linearity of the integral equation for 
the gluon Green's function. We will see below that Eq. ([5]) holds very similarly 
at NLO. This fact has been useful in the study of different NLO BFKL cross 
sections using numerical techniques in recent years (see Ref. [6]). 

After this brief introduction to the structure of BFKL cross sections and 
its iterative expression we now turn to the NLO case. The factorization into 
impact factors and Green's function will remain, while the kernel and trajectory 
will be more complex than at LO. We discuss these points in the next section. 



3 Different contributions at NLO 

To discuss the various contributions to NLO BFKL cross sections we follow 
Ref. [5]. We comment in more detail those points which will turn out to be 
more relevant for our later discussion of inclusive jet production. Our starting 
point are Eqs. (jT|) to ((4]), which remain unchanged. Since at NLO the sr scale 
is no longer a free parameter, we should modify Eq. ([3J to read 



•^AB+ll __ r (sH;0,l) 



2 s 



1 A 



i /jn^r" +i r ^,„ +1 ), (9) 

q n+ l \SR;, h n+lJ 

The propagation of a Reggeized gluon with momentum qi in MRK takes place 
between two emissions with momenta ki-i and ki (see Fig. [!}. Therefore, at 
NLO, the term sr, which scales the invariant energy Sj-i.i, does depend on 
these two consecutive emissions and, in general, will be written as s_R ; i-i,i. It 
is important to note that the production amplitudes should be independent of 
the energy scale chosen and, therefore, 
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S R;n,n+l 



for the particle-particle-Reggeon vertices and 
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for the Reggeon-Reggeon-gluon production vertices. 

At NLO, besides the two-loop corrections to the gluon Regge trajectory, 
there are four other contributions which affect the real emission vertex. The 
first one consists of virtual corrections to the one gluon production vertex. The 
second stems from the fact that in a chain of emissions widely separated in 
rapidity two of them are allowed to be nearby in this variable, this is known 
as quasi-multi- Regge kinematics (QMRK). A third source is obtained by per- 
turbatively expanding the Reggeon propagators in Eq. §§§ while keeping MRK 
and every vertex at LO. A final fourth contribution is that of the production of 
quark-antiquark pairs. The common feature of all of these new NLO elements 
is that they generate an extra power in the coupling constant without building 
up a corresponding logarithm of energy so that a s (a s lns) n terms are taken 
into account. 

With the idea of introducing a jet definition later on, it is important to 
understand the properties of the production vertex which we now describe in 
some detail. 

Lets start with the virtual corrections to the single-gluon emission vertex. 
These are rather simple and correspond to Eq. f8} with the insertion of a single 
kernel or impact factor with NLO virtual contributions (noted as (v)) while 
leaving the rest of the expression at Born level (written as (£?)). More explicitly: 
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n4 fl) (q*.*+i) 

i=l 



II 4 B) (qi,qi+i; 

i=j+l 




(12) 



Now we turn to the discussion of how to define QMRK. For this purpose 
the introduction of an extra scale is mandatory in order to define a separation 
in rapidity space between different emissions. As in Ref. [5] we call this new 
scale s\- At LO MRK implies that all s,j = (fcj + k 3 ) 2 are larger than s\. In 
rapidity space this means that their rapidity difference \yi — yj\ is larger than 
In fj, . As we stated earlier, in QMRK one single pair of emissions is allowed 



to be close in rapidity. When any of these two emissions is one of the external 
particles A or B it contributes as a real correction to the corresponding impact 
factor. If this is not the case it qualifies as a real correction to the kernel. This 
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is summarized in the following expression where we write real corrections to the 
impact factors as (r): 
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n4 B )( qi , qi+1 ) 
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JJ 4 B )(qi,q i+1 
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a 2 



^-qmrk (qj i q^'+i ) 

(B) , 

(q 



n+l, 
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(13) 



The modifications due to QMRK belonging to the kernel or to the impact factors 



are, respectively, /Cqmrk and $p , i 



A^QMRK(q 4 , qi+i) 
4 r) (k) 



c/.s 



1 / ds 



Irr vrr-+gg(s) 0(sa - s) 

(2 7 r)°q 2 q 2 +1 
^p/?. vpr->pg(s) 0(sa - s) 



(14) 



(15) 



(2tt)s 

In both cases s denotes the invariant mass of the two emissions in QMRK. The 
Heaviside functions are used to separate the regions of phase space where the 
emissions are at a relative rapidity separation smaller than sa- It is within this 
region where the LO emission kernel is modified. (Jrr^gg and <jpr—>pg are the 
total cross sections for two Reggeons into two gluons, and an external particle 
and a Reggeon into an external particle and a gluon, respectively. / stands for 
the invariant flux and N c for the number of colours. 

For those sectors remaining in the MRK we use a Heaviside function to 
keep Si-i.i > sa, in this way MRK is clearly separated from QMRK. We then 
follow the same steps as at LO and use Eq. ([7]) with the modifications already 
introduced in Eq. ©, i.e. 

-(n+l) , n 
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After performing the integration over the Sj_i^ variables the following interest- 
ing dependence on s\ arises: 



(n+l) 
MRK 
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n+l 



n^ B )( qi ,q l+1 ) 
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4^+2) 

q«+2 



(17) 



It is now convenient to go back to Eq. ([I]) and write the lower limit sq of 
the Mellin transform as a generic product of two scales related to the external 
impact factors, i.e. sq = ^/so-a sq-b- By expanding in a s the factors with powers 
in uj and u>i it is then possible to identify the NLO terms: 
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(18) 



To combine this expression with that of the QMRK contribution we should 
make a choice for sr. The most convenient one is sr.jj = *Js~h~i s R-,ji where 
for intermediate Reggeon propagation we use sr-i = k 2 , and for the connection 
with the external particles sr-q = sq-a and s.R ; n+2 = sq-b- We can then write 



(n+l) 



MRK 



n 



(a; — 2u;i) , 



[± (w - 2uji) 
In 



4 5) (qi) 
q? 



n+l 



n4 B) (q,q. + i) 



,»=i 



q«+2 



klS ;A 



n+l r 

E 



(w - 2cjj) s 
in 



(w - 2ui n+2 ) 



ktik? 



In 



k n+l S 0;B 



(19) 



This corresponds to the LO result for J^( n+1 ) plus additional terms where the 
uj — 2u>i factor cancels, in such a way that they can be combined with the LO 
result of T {n \ 

The quark contribution can be included in a straightforward manner since 
between the quark-antiquark emissions there is no propagation of a Reggeized 
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gluon. In this way one can simply write 
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(20) 



The production kernel can be written as 

/C Q0 (q,,q 4+1 ) = {N 2 C ~1) I ds 



IRR a RR-^QQ(s) 

(2^q?q? +1 ' 



(21) 



with ct rr _^qq being the total cross section for two Reggeons producing the 
quark-antiquark pair with an invariant mass s. 

The combination of all the NLO contributions together generates the follow- 
ing expression for the NLO cross section: 
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(22) 



where the NLO real emission kernel contains several terms: 
Mq 4 ,q t+1 ) = f/C^+ZC^Wq.+r) 



= (}C^+]C^+}C GG +}C Q Q)(q i ,q i+1 ), (23) 



with ICqq given by Eq. (|21[) . The two gluon production kernel ICqg is the 
combination of /Cqmrk of Eq. (fT4")) and the MRK contribution in Eq. p^|) . It 
explicitly reads 



^GG(q»,q 4 +i) = {Nl - 1) / ds 



JrR<7RR->Gg(s) #( s A - S) 



(2^q 2 q 2 +1 



■/ d D ~ 2 q 4 B )( qi ,q)/C( S )(q,q l+1 )iln 



(q; - q) 2 (q»+i - q) 



(24) 



Below we will show that when s\ is taken to infinity the second term of this ex- 
pression subtracts the logarithmic divergence of the first one. When computing 
the total cross section it is natural to remove the dependence on the parameter 
sa in this way. For our jet production cross section, however, we prefer to retain 
the dependence upon sa. 
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For the impact factors a similar expression including virtual and MRK cor- 
rections as in Eq. (fT5|) arises: 



/ ^ 2 q $(f)(q)/C( s )(q, qi )iln ( ? i^—) . (25) 

j 2 v (qi - q) s o ; p / 

From this expression it is now clear why to choose the factorized form sq = 
y/SQ-A sq-b: in this way each of the impact factors <&a,b carry its own sq_a.b 
term at NLO independently of the choice of scale in the other. 

To conclude this section, for the sake of clarity, the different contributions 
to the NLO BFKL kernel 
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are pictorially represented in Fig. [2] 




(a) (5) (c) (d) (e) 



Figure 2: Contributions to real emission kernel at LO (a) and NLO (b-e). 



As a final remark we would like to indicate that the divergences present in the 
gluon trajectories (see Ref. [2]) are all cancelled inside the inclusive terms. 
We will see how the soft and collinear divergences of the production vertex are 
either cancelled amongst its different components or are regularized by the jet 
definition. 

After having introduced the notation and highlighted the different con- 
stituents of a BFKL production kernel at NLO, in the coming section we describe 
how to calculate the inclusive production of jets in two different environments. 
The first one is the case of the interaction between two small and perturba- 
tive objects, highly virtual photons, and the second will be the collision of two 
large and non-perturbative external particles such as the ones taking place at 
hadron-hadron colliders. 
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4 Inclusive jet production at LO 



As MRK relies on the transverse scales of the emissions and internal lines being 
of the same order it is natural to think that processes characterized by two large 
and similar transverse momenta are the ideal environment for BFKL dynamics 
to show up. Moreover, as the resummation is based on perturbative degrees of 
freedom, these large scales associated to the external particles should favor the 
accuracy of the predictions. An ideal scenario is the interaction between two 
photons with large virtualities Q\ 2 in the Regge limit s ^> \t\ ~ Q\ ~ Q\. The 
total cross section for this process has been investigated in a large number of 
papers in recent years. Here we are interested in the inclusive production of 
a single jet in the central region of rapidity in this process. We will consider 
the case where the transverse momentum of the jet is of the same order as the 
virtualities of the photons. 

As a starting point we review single jet production at LO accuracy. As 
usual the total cross section can be written as a convolution of the photon 
impact factors with the gluon Green's function, i.e. 

f d 2 k a f d 2 k b ^ n , ^ „ , f s+toc du ( sY „ n , , ,„„, 
J 2irk 2 a J 2irk 2 J g _ ioo 2m \s J 

A common choice for the energy scale is Sq = |k a | |kt| which naturally introduces 
the rapidities y i and of the emitted particles with momenta and since 

^— j = e "(vA-va). (27) 

Let us remark that a change in this scale can be treated as a redefinition of 
the impact factors and, if sq is chosen to depend only on k Q or only on kb, the 
kernel as well. This treatment lies beyond LO and will be discussed in the next 
section. The gluon Green's function f u corresponds to the solution of the BFKL 
equation 

iufjk a ,k b ) = 5 2 (k a -k fc )+ J d 2 k/C(k Q ,k)/ w (k,k b ), (28) 

/C(k a ,k) = 2^(k 2 )5 2 (k a -k)+/C r (k a ,k), (29) 

where the kernel K. contains a term related to the Reggeized gluon propagator, 
the trajectory w(k 2 ), and the real emission kernel, K r . 

For the inclusive production of a single jet we assign to it a rapidity yj 
and a transverse momentum kj, as shown in Fig. [3J In this way, if kj = 
ctjPA + PjPb + kj± the corresponding rapidity is yj — | In Using its on- 
shell condition we can write 

kj = \^-e yj p A + J^e-y j p B + kj±. (30) 
V s V s 
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Figure 3: Total cross section and inclusive one jet production in the BFKL 
approach. 



It is possible to single out one gluon emission by extracting its emission 
probability from the BFKL kernel. The differential cross section in terms of the 
jet variables can then be constructed in the following way: 

/ 535/535 •**■>••<*> 

'/MMr^)^--' 

v/ 1 ^ f 5+100 dcu' fs BJ y' 

X V(q ,q b ;kj,yj) / — — L>(-<lb,-k b ) (31) 

with the LO emission vertex being 

V(q 0) q 6 ; k J; yj) = K,f^ (q a , -q 6 ) 5^ (q a + q b - kj) . (32) 

By selecting one emission to be exclusive we have factorized the gluon Green's 
function into two components. Each of them connects one of the external par- 
ticles to the jet vertex. In the notation of Eq. (|31|) the energies of these blocks 
are 

sa.j =(pa + qb) 2 , s B .j ={pb + q a ) 2 ■ (33) 

In a symmetric situation, where the jet provides a hard scale as well as the 
impact factors, a natural choice for the scales is similar to that in the total cross 



da 
d 2 kjdyj 
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section 



|k„||k. 



So=|kj||k 6 



(34) 



These choices can now be related to the relative rapidity between the jet and 
the external particles. To set the ground for the NLO discussion of the next 
section we introduce an additional integration over the rapidity rj of the central 
system: 



da 



(Pkjdyj 



d 2 k a 

27Tk2 



d q a / d q b / di] 



du> 
2tH' 



$A(k a ) 



(5+zoo 
5—ioo 
b 



,u(yA-v) 



8+ioo 



/w(k a , qa) 



V(qa,q ,??;k,7,yj) 



ifcrkg 



/^'(-q6,-k b ) (35) 



with the LO emission vertex being 

V(q a , q 6 , m kj, yj) = /C< B > (q a , -q 6 ) (q a + q 6 - kj) 6( V - y,j). (36) 

Eqs. (|35[) and (f3"6")l will be the starting point for the NLO jet production in the 
symmetric configurations. 

Let us now switch to the asymmetric case. In general we can write q a and 
Qb as 



q a =a a PA + PaPB + q a ± 



q b =a b p A + (3 b p B + qb±- 



(37) 



The strong ordering in the rapidity of emissions translates into the conditions 
a a ^> a and (3b ^> I3 a . This, together with momentum conservation q a + q b = 
kj, leads us to aj — a a + a b ~ a a , flj — Pa + 0b ~ Pb and 



SAJ =PjS, 



sbj =ajs. 



(38) 



While the longitudinal momentum of q a {qb) is a linear combination of pa and 
Pb we see that only its component along pa(pb) matters. 

If the colliding external particles provide no perturbative scale as it is the 
case in hadron-hadron collisions, then the jet is the only hard scale in the process 
and we have to deal with an asymmetric situation. Thus the scales sq and s 
should be chosen as k % alone. At LO accuracy sq is arbitrary and we are indeed 
free to make this choice. Then the arguments of the gluon Green's functions 
can be written as 



SAJ 

so 



1 

Ola 



SBJ 

so 



1 

'Tb 



(39) 



The description in terms of these longitudinal components is particularly useful 
if one is interested in jet production in a hadronic environment. Here one can 
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introduce the concept of unintegrated gluon density in the hadron. This rep- 
resents the probability of resolving a gluon carrying a longitudinal momentum 
fraction x from the incoming hadron, and with a certain transverse momentum 
fcx- With the help of Eq. ([55)1 a LO unintegrated gluon distribution g can be 
defined from Eq. (p?Tj) as 



g(x,k) 



d\ 



2vrq 2 

Then we can rewrite Eq. (|31[) as 
da 



$ P (q) 



5+ioo 



S-i 



2ttI 



^/ u (q,k). 



d q Q / efei d q b dx 2 



x q a ).g(a;2, q&)V(q a , an, q&, x 2 ; kj, yj) 
with the LO jet vertex for the asymmetric situation being 

V(q, a ,xi,q b ,x 2 ;kj,yj) = K.^ (q Q ,-q b ) 



5 (2) (q Q + q b - kj) <5 L - y^e^j <5 ^ 2 - 




(40) 



(41) 



(42) 



Having presented our framework for the LO case, in both 7*7* and hadron- 
hadron collisions, we now proceed to explain in detail what corrections are 
needed to define our cross sections at NLO. Special attention should be put into 
the treatment of those scales with do not enter the LO discussion but are crucial 
at higher orders. 



5 Inclusive jet production at NLO 

A similar approach to that shown in Section 4 remains valid when jet production 
is considered at NLO. The crucial step in this direction is to modify the LO jet 
vertex of Eq. (|36|) and Eq. (j42| to include new configurations present at NLO. We 
show how this is done in the following first subsection. In the second subsection 
we implement this vertex in the symmetric 7*7* case, and we repeat the steps 
from Eq. (|2l)|) to Eq. ([55]) , carefully describing the choice of energy scale at each 
of the subchannels. In the third subsection hadron-hadron scattering is taken 
into consideration, and we extend the concept of unintegrated gluon density of 
Eq. (|4"D|) to NLO accuracy. Most importantly, it is shown that a correct choice 
of intermediate energy scales in this case implies a modification of the impact 
factors, the jet vertex, and the evolution kernel. 



5.1 The NLO jet vertex 

For those parts of the NLO kernel responsible for one gluon production we 
proceed in exactly the same way as at LO. The treatment of those terms related 
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to two particle production is more complicated since for them it is necessary 
to introduce a jet algorithm. In general terms, if the two emissions generated 
by the kernel are nearby in phase space they will be considered as one single 
jet, otherwise one of them will be identified as the jet whereas the other will 
be absorbed as an untagged inclusive contribution. Hadronization effects in 
the final state are neglected and we simply define a cone of radius Rq in the 
rapidity-azimuthal angle space such that two particles form a single jet if i?i 2 = 
yj (</>i — </> 2 ) 2 + (yi — y-i) 2 < Rq- As long as only two emissions are involved this 
is equivalent to the fcr-clustering algorithm. 

To introduce the jet definition in the 2 — > 2 components of the kernel it is 
convenient to start by considering the gluon and quark matrix elements together: 



(/Cqmrk +/C Q q) (q a , -q 6 ) = J d D 2 k 2 J 



dy 



2 



x ( |^ 2 g(qa,q6,ki,k 2 )| 2 9(sa - S12) + |A g (qa,qf,,ki,k 2 )| 2 J , (43) 

with A2P being the two particle production amplitudes of which only the gluonic 
one also contributes to MRK. This is why a step function is needed to separate 
it from MRK. Momentum conservation implies that k x = q a + q^ — k 2 . 

The expression (|43[) is not complete as it stands since we should also include 
the MRK contribution as it was previously done in Eq. (f2"4")) : 

(K, GG + K, Q q) (q , -q 6 ) = J d D ~ 2 k 2 J dy 2 \B(q a , q 6 , ki, k 2 )| 2 
= J d D - 2 k 2 J d2/ 2 ||A 9 (qa,q6,ki,k 2 )| 2 e(s A -si 2 ) 

- (q Q , q Q - ki) /C< s )(q a - k x , -q 6 ) ±6 (hx^ - y^j [y 2 - In ^) 
+ |^l 29 (q Q ,q b ,k 1 ,k 2 )| 2 |. (44) 
We are now ready to introduce the jet definition for the double emissions. The 



Ki 



NLO versions of Eq. ([55)) and Eq. (|42[) then read, respectively, 



[v] 

(a) 



V(q a , q 6 , m kj, yj) = ^ + fc< >) (q a , -q 6 ) 

+ J d°- 2 k 2 dy 2 |£(q Q ,q b ,k, 7 -k 2 ,k 2 )| 2 #(i? - i? 12 
+ 2 /" ^- 2 k 2 dyalBfqa.qj.kj.ka)! 2 9(R J2 - Rq) 
V(q a ,x 1 ,q b ,x 2 ]'kj,yj) = (IC { r B) + ) (qo, -Q6 



(45) 



+ / d u ~ 2 k 2 dy 2 \B(q a ,q b ,kj -k 2 ,k 2 )| Z 6(R - R12, 



\x] 
(c)' 



(>') 



+ 2 / d u ~ z k 2 dy 2 \B{q a ,q b ,kjM)\ 0{R J2 ~ R ) 

In these two expressions we have introduced the notation 

[y] 



(a,b) 

ly] 

(c) 

[«] 

(o,6) 
[x] 

(c) 



<5 (2) (q a +q 6 -k, 7 )%-y^ b )), 



<5 (2) (q a +q 6 -k 7 -k 2 )(5(?7-j/ 



(o,6) 



<5 (2) (q»+q6-kj)*(xi-a!i ' 6) )*(x2 
<5 (2) (q a + - kj - k 2 ) <5 (xi - x[ c) ^j S (x 2 - x. 



(46) 

(47) 
(48) 
(49) 
(50) 



The various jet configurations demand different y and x configurations. These 
are related to the properties of the produced jet in different ways depending on 
the origin of the jet: if only one gluon was produced in MRK this corresponds 
to the configuration (a) in the table below, if two particles in QMRK form a jet 
then we have the case (b), and finally case (c) if the jet is produced out of one of 
the partons in QMRK. The factor of 2 in the last term of Eq. (gSJ) and Eq. (|4l)j) 
accounts for the possibility that either emitted particle can form the jet. Just 
by kinematics we get the explicit expressions for the different x configurations 
listed in the following table: 



JET y configurations 



x configurations 



a, 
b) 
c) 



y [a) = yj 
y {h) = yj 



r (a) _ \kj\ e y., 
1 



(a) _ Jkj| yj 
Jb) _ y/S yj 

X 2 - v /i e 



(c) _ \kj[ e yj , \k2[ e y 2 X (C) _ 



hjl^-yj _|_ k 2l r -vi 



oft 



0T" 



The variable £ is defined below in Eq. (|58|) . Due to the analogue treatment 
of the emission vertex either expressed in terms of rapidities or longitudinal 
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momentum fractions in the remaining of this section we will imply the same 
analysis for both. In particular, we will not explicitly mention these arguments 
when we come to Eqs. ([551 loU)) . 

The introduction of the jet definition divides the phase space into different 
sectors. It is now needed to show that the final result is indeed free of any 
infrared divergences. In the following we proceed to independently calculate 
several contributions to the kernel to be able, in this way, to study its singularity 
structure. 

The NLO virtual correction to the one-gluon emission kernel, fc( v \ was 
originally calculated in Ref. [7J [5J [H] • Its expression reads 



IC { r v) (q a ,-q 6 ) 
N r . e 



T(l - e) A 2 



7T- 1 



q 2 . 



3A 2 /q 2 
2 m ~2 

q 2 . Wb 



A 2 

M 2 
si 

q 2 . 



i 



7T" 

T 



2eC(3) 



1 



Nr 



q 2 



4 



-In 

A 2 (q 2 +q^ -4q a q 6 ) 



3(q 2 "^) 



2\2 



6qaq& 



iq a - qtj 



A 4 (q 2 



■<£) 



6qaqb(q 2 - q& 



2V2 



(qa + q6-4q Q qb) 




(51) 



with /3q = (HN C — 2n/)/3, ((n) — YlkLi ^ ™ an( I A = q a + q^. can be 
expressed in terms of the renormalized coupling constant g^ in the MS renor- 
malization scheme by the relation g 2 = <? 2 A r c T(l — e) (47r)~ 2 ~ e . Note that the 
expression for the virtual contribution given in [4] lacks the log squared. 

Those pieces related to two-gluon production in QMRK can be rewritten in 
terms of their corresponding matrix elements as 

^QMRK(q a , -qt) = J d D ~ 2 k 2 J dy 2 |„4 2g (q a , q&, k i; k 2 )| 2 6>(s A - s 12 ) 

<? 2 /i" 2e iV 2 f d D ~ 2 k 2 f 
= ^p^f J ^(2ttP-4 J dV2 Aglu ° ns 9(SA ~ Sl2) ' (52) 

and those related to quark-antiquark production are 



^oo(qa, -q&) 



2 k 2 / dy 2 |A g (qa,qf,,ki,k 2 )|' 



gl^Nl 



dy 2 A 



quarks ■ 



(53) 



^(2^) c + 1 q 2 q 2 7 /^(2tt) c - 4 
We have calculated the corresponding amplitudes, using the Mandelstam invari 
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ants s, t, and u, and our results are 

,2„2 



n 2 n 2l 1 , 1 « 1 I"* 1 

q a q& S ~ 77 



iu itu kj'k 2 . 4 V cc k|t 1 - x k 2 u / 4k 2 k2. 











♦G 


4) 








1 - 


x 1 




X 


k 2 

K 2 





1 - X 



q 2 . 



1-ac 



4 V s 



1 



si 

4s 



3- k 2 , 

i - xicf 



1 - xk( 



2 1 



D-2j /(kx-qj^^-q.) 2 -^ 



4 | \ it! 
1 / fk. - q Q ) 2 - ^k 2 £ A ^ (k] _ qa ) 2 _ l^ k 2 ^ ^ 

t s 



.,(54) 



iquarks 



,2_2 











(l- 


1) 









: k? 



l-a; 



- q ) 2 (k 2 - q a ) 2 - kfkf 2 



f u 



1 /(k 2 -q a ) 2 



g l f 2 



E\ /(k!-q a ) 2 -^k 2 £ 



47V? 



(ki-q^^-q^-kfk 2 ^ - 2 ~ 2 



tu 



(55) 



These expressions are in agreement with the corresponding ones obtained in 
Ref. [1]. The following notation has been used: 



|ki| 



A 

E 

E 



|ki| + |k 2 |e^' 
(1 - x)ki - xk 2 , 
A 2 



s + A 



x(l — x) 



1 — 2x 

2(2x - l)q 2 + 4Aq a + — T A 2 

x(l — X) 

2x(l - x) ((2x - 1) A 2 + 2AA) 



q 2 , 



x(l - x)A 2 + A 2 ' 



(56) 
(57) 
(58) 

(59) 



We now study those terms which contribute to generate soft and collinear di- 
vergences after integration over the two-particle phase space. They should be 
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able to cancel the e poles of the virtual contributions in Eq. (f5Tj) , i.e. 



Here we identify those pieces responsible for the generation of these poles. 

One of the divergent regions is defined by the two emissions with momenta 
ki = aipA + PiPb + k\± and fc 2 = ctiVA + fiiVB + &2_l becoming collinear. This 
means that, for a real parameter A, k\ ~ A&2, fcij_ — Afc2±, ol\ ~ A«2 and 
thus o^fcii. — «ifc2J_ — 0. Since cc = a this is equivalent to the condition 
A ~ 0. In the collinear region s = x (^_ x ^ tends to zero and the dominant 
contributions which are purely collinear are 



^singula 



* singula] 
quarks 



= _ q|q|2 

collinear S S 

collinear 2N C sE 



D-2E 2 
16 I 2 " 

8N C s 2 ' 



A { i)+A 



(2)' 



(61) 
(62) 



The quark-antiquark production does not generate divergences when ki or k2 
become soft, therefore we have that the only purely soft divergence is 



.4 



singular 
;luons 



2 2 i 



1 



2 2 



1 



Atu kfk 2 . 



4k 2 k 2 



.4 



(3) 



.4 



(4) 



2A 



(4). 



(63) 



where we have used the property that, in the soft limit, the tu product tends 
to q 2 q 2 . We will see that these terms will be responsible for simple poles in e. 
The double poles will be generated by the regions with simultaneous soft and 
collinear divergences. They are only present in the gluon-gluon production case 
and can be written as 



A 



singular 



2 2 

lilt 

4s£ 



soft&collincar 

l 



2 2 

4s 



x 1 



1 



k 2 

K 2 



-. — tt + — + q a 

(1 — x)t xu 



k 2 

K 2 



xk\ 
1 



1 ~xk( 



t 



A 



(5) 



.4 



(6)- 



l-xkj 



k 2 

K 2 



(64) 



Focusing on the divergent structure it turns out that in the soft and collinear 
region the first line of Eq. (ffi4")> , , has exactly the same limit as the second 
line, Arty . This is very convenient since we can then simply write 



.4 



singular 



2 2 



1 - X 1 



soft&collincar 2s V X ko 1 

The MRK contribution of Eq. (gH) has the form A M rk 



xk\ 



2A 



(5)- 



(65) 



-4A(4) and when 
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added to all the other singular terms we get the expression 
J d°- 2 k 2 J dy 2 \B s (q a ,q b ,k 2 M)\ 2 = 



with 

0(s A - Sl2 ) + A» 



dV2 {^3„'r 0(sa - Sl2 ) + (66) 



q 2 q 2 2 | D-2E' 
S 5 16 s 2 



2 2 

2k 1 k 2 



2.1 



k 2 

K 2 



+ 



Gluon| coUl GluonUu 2 Gluon| aoft 



Gluon| Bot t&coll 



nj E 2 




0{SA - 8l 2 ) 



(67) 



We have labeled the different terms to study how each of them produces the e 
poles. We will do this in Section 6. 

With the singularity structure well identified we now come back to Eqs. (l45l 
46|) and show how they are free of any divergences. Only if the divergent terms 
belong to the same configuration this cancellation can be shown analytically. 
With this in mind we add the singular parts of the two particle production of 
Eq. (|6"6"j) in the configuration (a) multiplied by = 1 — 8(Ro — R12) — 9{Ri2 — Ro)' 



V 



K. 



d D - 2 k 2 dy 2 



J d D 2 k 2 dy 2 |6 s (q a ,q 6 ,kj - k 2 ,k 2 



(a) 



|i3(q a ,q b ,k,7 - k 2 ,k 2 ) 



(!>) 



(a) 



- |S s (qa,qb,kj - k 2 ,k 2 )| 
2 J d°- 2 k 2 dy 2 |£(q a ,q b ,k 7 ,k 2 )| 2 0(.R, 72 -i? o 

- J d°- 2 k 2 dy 2 \B s (q a , q b , k, - k 2 , k 2 )| 2 0{R 12 - R ) 



9(Ro — R12) 



(a) 



(68) 



The cancellation of divergences within the first two lines is now the same as 
in the calculation of the full NLO kernel. In Section 6 we will show how the 
first two lines of Eq. (|6"5|) are free of any singularities in the form of e poles. In 
doing so we will go into the details of the role of sa- The third and fourth lines 
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are also explicitly free of divergences since these have been subtracted out. The 
sixth line has a ki <-> k 2 symmetry which allows us to write 



V = 



(B) 



d D ~ 2 k 2 dy 2 



•/CW)(q a) - qi) ) 



J d D 2 k 2 dy 2 |£ s (q Q ,q b ,kj - k 2 ,k 2 ) 



(a) 



|£(q a ,q b ,kj - k 2 ,k 2 )|' ! 



(i>) 



|£ s (q a ,q b ,kj - k 2) k 2 )| i 



(a) 



0(Ro — R12) 



2 / d D ~ 2 k 2 dy 2 



\B(q a ,q b ,kj,k 2 )\ 2 6(R J2 -R Q : 



|6 s (q a ,q 6 ,k, / -k 2 ,k 2 )| 2 0(i? 12 -i? o )0(|k 1 | - |k 2 | 



(a) 



(69) 



We can now see that the remaining possible divergent regions of the last line 
are regulated by the cone radius i?o- 

It is worth noting that, apart from an overall a 2 (/i 2 ) factor, the NLO terms in 
the last four lines in Eq. (|69p do not carry any renormalization scale dependence 
since they are finite when e is set to zero. The situation is different for the first 
two lines since V contains a logarithm of /i 2 in the form 



V = V< B > ( 1 



a s{y 2 ) §0 
Air N r 



k 2 

A 1 



AV. 



(70) 



where AV contains the third to sixth lines and the /i-independent part of the 
first two lines of (1691). It is then natural to absorb this term in a redefinition 



of the running of the coupling and replace a s (/j, 2 ) by a s (k 2 ). For a explicit 
derivation of this term we refer the reader to Section 6. 

Therefore we now have a finite expression for the jet vertex suitable for 
numerical integration. This numerical analysis will be performed elsewhere since 
here we are mainly concerned with the formal introduction of the jet definition 
and the correct separation of the different contributions to the kernel. 

What remains to be proven is the cancellation of divergences between Eq. ([60]) 
and Eq. (|66|) . This will be performed in Section 6. Before doing so, in the next 
two subsections, we indicate how to introduce our vertex in the definition of the 
differential cross section. Especial care must be taken in the treatment of the 
energy scale in the Reggeized gluon propagators since in the symmetric case it 
is directly related to the rapidity difference between subsequent emissions, as we 
will show in the next subsection, but in the asymmetric case of hadron-hadron 
collisions it depends on the longitudinal momentum fractions of the t-channel 
Reggeons. 
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5.2 Production of jets in 7*7* scattering 



We now have all the ingredients required to describe the inclusive single jet 
production in a symmetric process at NLO. To be definite, we consider 7*7* 
scattering with the virtualities of the two photons being large and of the same 
order. All we need is to take Eq. (|3"Tj) for the differential cross section as a 
function of the transverse momentum and rapidity of the jet. The vertex V to 
be used is that of Eq. in the representation based on rapidity variables of 
Eq. (I45p . The rapidities of the emitted particles are the natural variables to 
characterize the partonic evolution and s-channel production since we assume 
that all transverse momenta are of the same order. 

Let us note that the rapidity difference between two emissions can be written 

as 

V% - Ui+i = In ' S " +1 (71) 
\/ k i k i+i 

which supports the choice sr^^+i = ^/k|k| +1 in Eq. This is also technically 
more convenient since it simplifies the final expression for the cross section in 

Eq. (m. 

In Fig. [4] we illustrate the different scales participating in the scattering 
and the variables of evolution. We write down the conditions for MRK: all 
transverse momenta are of similar size and much larger than the confining scale, 
the rapidities are strongly ordered in the evolution from one external particle 
to the other. At each stage of the evolution the propagation of the Reggeized 
gluons, which generates rapidity gaps, takes place between two real emissions. 
There are many configurations contributing to the differential cross section, each 
of them with a different weight. Eq. pip represents the sum of these production 
processes. 



5.3 The unintegrated gluon density and jet production in 
hadron— hadron collisions 

In this subsection we now turn to the case of hadron collisions where MRK has 
to be necessarily modified to include some evolution in the transverse momenta, 
since the momentum of the jet will be much larger than the typical transverse 
scale associated to the hadron. 

In the LO case we have already explained that, in order to move from the 
symmetric case to the asymmetric one, it is needed to change the energy scale 
from Eq. (|34[) to Eq. (j3"5|) . This is equivalent to changing the description of 
the evolution in terms of rapidity differences between emissions to longitudinal 
momentum fractions of the Reggeized gluons in the i-channel. Whereas in 
LO this change of scales has no consequences, in NLO accuracy it leads to 
modifications, not only of the jet emission vertex but also of the evolution 



23 



MRK: p\ ~ p| ~ k 2 ~ k 2 > A2j CD , ^ = j/o > j/i > ■ ■ • > y n > t/ n+1 
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Figure 4: Momenta for 2 ^ 2 + (n — 1) + jet amplitude in the symmetric 
configuration with MRK. The produced jet has rapidity yj = yj and transverse 



momentum kj = k. 
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Figure 5: Momenta for 2 — > 2 + (n — 1) + jet amplitude in the asymmetric 
configuration with fc t -ordered MRK. 
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kernels above and below the jet vertex. These new definitions will allow the 
cross section still to be written in a factorizable way and the evolution of the 
gluon Green's function still to be described by an integral equation. 

To understand this in detail we start by writing the solution to the NLO 
BFKL equation iteratively, i.e. 



I 



d 2 k a / w (k OJ q a ) 



OO 

E 



j ~ 1 f 1 



(72) 



where qi = k a and q^ = q a . We now focus on one side of the evolution towards 
the hard scale since the other side is similar and use Fig. |3] as a graphical 
reference. Starting with the symmetric case the differential cross section for jet 
production contains the following evolution between particle A and the jet: 



da 



d 2 kjdy,j 



d 2 q a / d 2 k a 



<Mko) 
27rk2 



/ 2^^ ( ' ka ' qa ' ) ( "Tlfp?) V ( qa ' qb;kj ' yj )' 



(73) 



In the asymmetric situation where k 2 3> k 2 the scale i/k^kj should be 
replaced by k 2 . In order to do so we rewrite the term related to the choice of 
energy scale. To be consistent with Fig. [5] we take kj = kj, ko = — k a = qi 
and qj = q Q . To start with it is convenient to introduce a chain of scale changes 
in every kernel: 



SAJ 



n 



,2 \ 2 



k 2 



SAJ 
k 2 7 



which can alternatively be written in terms of the i-channel momenta as 

n?-i 



(74) 
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SAJ 
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(75) 



For completeness note that we are indeed changing the variable of evolution 
from a difference in rapidity: 



SAJ 



va-vj 



(76) 



to the inverse of the longitudinal momentum fraction, i. e. 



SAJ 

k 2 



1 

a J 



(77) 



This shift in scales translates into the following expression for the cross 
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section: 
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As we mentioned above these changes can be absorbed at NLO in the kernels 
and impact factors, we just need to perturbatively expand the integrand. The 
impact factors get one single contribution, as can be seen in Fig. and they 
explicitly change as 



<f(k Q ) = $(k a )-ik2 



, 2q ^K(s) (q)ka) i n s! 

q 2 k2 



(79) 



The kernels in the evolution receive a double contribution from the different 
energy scale choices of both the incoming and outgoing Reggeons (see Fig. [5]). 
This amounts to the following correction: 

K( qi ,q 2 ) = /C( qi ,q 2 )-i /d 2 q/C( B )( qi ,q)/C( B )(q,q 2 )ln^. (80) 

* j q 2 

There is a different type of term in the case of the emission vertex where the 
jet is defined. This correction has also two contributions originated at the two 
different evolution chains from the hadrons A and B. Its expression is 

V(q a ,q b ) = V(q a ,q b )-iyrf 2 q/C( S )(q Q ,q)V( S )(q,q b )ln — 



q 6 ) 2 



\ f d 2 qV (B) (q a ,q)/C^(q,q b )ln ^ 

2 7 (q a - q) 



(81) 



These are all the modifications we need to be able to write our differential 
cross section for the asymmetric case. The final expression is 
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As in the LO case, we can use Eq. ([77| to define the NLO unintegrated gluon 
density as 



?(x,k) 



f , 2 Mg) f dLo ~ 

J dq -2^J ftri /w(k,q)aJ 



(83) 



The gluon Green's function f u is the solution to a new BFKL equation with the 
modified kernel of Eq. ([80)) which includes the energy shift at NLO, i.e. 



(84) 
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In this way the unintegrated gluon distribution follows the evolution equation 



dg(x,g.g) 
dlnl/x 



d 2 q/C(q a ,q).g(a;,q). 



(85) 



Finally, taking into account the evolution from the other hadron, the differential 
cross section reads 



with the emission vertex taken from Eq. . 

We would like to indicate that with the prescription derived in this subsection 
we managed to express the new kernels, emission vertex and impact factors as 
functions of their incoming momenta only. It is also worth mentioning that 
the proton impact factor contains non-perturbative physics which can only be 
modeled by, e.g. 



where Pi are positive free parameters, with Qq representing a momentum scale 
of the order of the confinement scale. The initial x dependence in this expression 
would be of non-perturbative origin. 

Let us also point out that the prescription to modify the kernel as in Eq. (j80|) 
was originally suggested in the first paper of Ref. [5] in the context of deep 
inelastic scattering. This new kernel can be considered as the first term in an 
all orders perturbative expansion due to the change of scale. When all terms 
are included the kernel acquires improved convergence properties and matches 
collinear evolution. Details of this procedure can be found in Ref. [TT], where 
the collinear resummation was done in Mellin space. In a future publication we 
intend to investigate how these corrections can be phrased in momentum space, 
and how they affect the behaviour of the unintegrated gluon distribution. For 
this we will use the procedure developed in Ref. [12] where the resummation to 
all orders corresponding to the energy shift was proven to be equivalent to a 
Bessel function of the first kind with argument depending on the strong coupling 
and a double logarithm of the ratio of transverse scales. 



6 Cancellation of divergences and a closer look 
at the separation between MRK and QMRK 



During the calculation of a NLO BFKL cross section, both at a fully inclusive 
level and at a more exclusive one, there is a need to separate the contribu- 
tions from MRK and QMRK. In order to do so we have followed Ref. [5] and 
introduced the parameter s\ in Eq. (Q3J) and Eq. (fT5")) . In principle, at NLO 
accuracy, our final results should not depend on this extra scale. In fact, as we 




! 2 q / d 2 q b g(x a ,q a ) g(x b ,q b ) V{q a ,q b ;kj,y.j) 7 (86) 




(87) 
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have remarked earlier in our discussion of the total cross section (after Eq. 
we could have taken the limit sa — * oo: the logarithms of sa cancel, and the 
corrections to the finite pieces die away as ©(s^ 1 ). In the context of the in- 
clusive cross section, however, we prefer to treat sa as a physical parameter: it 
separates MRK from QMRK and, hence, cannot be arbitrarily large. We will 
therefore retain the dependence upon sa- in the remainder of this section we 
demonstrate that, in our inclusive cross section, all logarithmic terms cancel 
(analogous to Eq. (|24p). and we will then leave the study of the corrections 
of the order 0(s^ 1 ) for a numerical analysis. It will also be interesting to see 
how this dependence on sa could be related to the rapidity veto introduced in 

Ref. peg. 

Let us consider the sa dependent terms in Eq. ([66]) which are only present 
in the gluon piece: 
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where we have used the numbering (7, //, III, IV) corresponding to, respec- 
tively, (Gluon| colll , Gluon| co n 2 , Gluon| soft , Gluon| so ft& co ii) in Eq. ([67j). 

To calculate each of the Si terms we start by transforming the rapidity 
integral into an integral over x in the form J dAy = J x ^fZ x ) ■ We consider 
sa much larger than any of the typical transverse momenta. In the limit of 
large sa the theta function 9(s\ — s) amounts to the limits + O (s^ 2 ) and 
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for the x integral. 
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We are only interested in the logarithmic dependence on sa and hence we do 
not need to calculate O (s^ 1 ) or sa independent factors. 
The next sa contribution we calculate is Siy which reads 
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The part with logarithmic s\ dependence can be calculated analytically: 
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(91) 



It is then clear that this logarithmic sa contribution cancels against that of Sm 
in Eq. (89]). 

Let us proceed now to show that the contribution of Si is directly of O (s^ 1 ) 
and does not contribute with any logarithm of sa. In the relevant integral we 
introduce the change of variable k 2 — > A = (1 — x) A — k 2 and obtain 
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We do not write here the lengthier but similar expression which corresponds to 
Sn and also only contributes to O (s^ 1 ). 

With this we have shown that the sum of different terms in Eq. ([88]) is free of 
logarithmic dependences on sa proving, in this way, that the remaining O (s^ 1 ) 
corrections vanish at large values of sa- In particular, it is possible to take the 
sa — * oo limit in order to completely eliminate the dependence on this scale. 
This is convenient in the fully inclusive case where it is very useful to write a 
Mellin transform in the kx dependence of the NLO BFKL kernel. 

If we perform this sa — > oo limit then Sm and Siv can be put together and 
their sum is 
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Regarding Si the integration gives us 
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The contribution from 5// is more complicated and the relevant integral can be 
obtained in the following way: 
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We now need to integrate it over x to obtain: 
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This result gives the same poles in e as the result given in [4], but differs for 
the finite contribution. To obtain all the e poles we now also include the quark 
contributions present in Eq. ([6"6"j) . We denote them as 
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where the correspondence with Eq. (f6"7)) is (V, V/) — > (Quark] coiij, Quark| co u 2 ). 
Adding everything up the sum of all the terms reads 
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The final expression for Eq. (JSSJ) is then 
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When we combine this result with the singular terms of Eq. (|51[) then we ex- 
plicitly prove the cancellation of any singularity in our subtraction procedure 
to introduce the jet definition. The finite remainder reads 
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We have already discussed the logarithmic term due to the running of the cou- 
pling in Eq. (|70[) . The non-logarithmic part is similar to that present in other 
calculations involving soft gluon resummations [14j where terms of the form 



a s (1 + S a s ) 



(102) 



appear and offer the possibility to change from the MS renormalization scheme 
to the so-called gluon-bremsstrahlung (GB) scheme by shifting the position of 
the Landau pole, i.e. 



At^ exp S— 
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The factor S differs from ours in the ir 2 term: 
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The origin of this discrepancy lies in the fact that we used the simplest form 
of subtraction procedure. In the Appendix we suggest a different subtraction 
term which is more complicated in the sense that it substracts a larger portion 
of the matrix element in addition to the infrared divergent pieces. When this is 
done and we put together the divergent pieces of Eq. (f5Tj) and the second line 
of Eq. (|119[) then we recover the same S term. 



7 Conclusions 

In this paper we have extended the NLO BFKL calculations to derive a NLO 
jet production vertex in fcT _ factorization. Our procedure was to 'deconstruct' 
the NLO BFKL kernel to introduce a jet definition at NLO in a consistent 
way. After a careful study of the different energy scales and contributions to 
the kernel we were able to show the infrared finiteness of this jet vertex and its 
dependence on the scale s\, which separates MRK from QMRK. As the central 
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result of this paper, we have defined the jet production vertex (|6T))) in terms of 
longitudinal momentum fractions, and we have explicitly given the necessary 
subtraction, both at the matrix element level (f5Tj) as well as integrated over 
the corresponding phase space (|100[) . Our calculations also suggest that the 
natural scale for the running of the coupling at the jet vertex is the square of 
the transverse momentum of the jet ([70)) . We have shown how this vertex can 
be used in the context of 7*7* or hadron-hadron scattering (f86| to calculate 
inclusive single jet cross sections. For this purpose we have formulated, on the 
basis of the NLO BFKL equation, a NLO unintegrated gluon density valid in 
the small-x regime. 

In our analysis we have been careful to retain the dependence upon the en- 
ergy scale sa which appears at NLO accuracy and separates multi-Regge kine- 
matics from quasi-multi-Regge kinematics. In the NLO calculation of the total 
cross section, one may be tempted to take the limit sa — > 00, thus disregarding 
the 1/sa corrections to the NLO BFKL kernel. However, when discussing inclu- 
sive (multi-) jet production one has to remember that sa has a concrete physical 
meaning: it denotes the lower cutoff of rapidity gaps and thus directly enters 
the rapidity distribution of multi-jet final states. In a self-consistent description 
then also the evolution of the unintegrated gluon density has to depend upon 
this scale. 

Hence we are well prepared for our next step, the numerical study of single 
or multiple jet production in hadron-hadron collisions at the LHC. One issue 
to be covered will be the question of handling the running of coupling. Further 
applications of our NLO fey-formalism include W and Z as well as heavy flavor 
production in the small-x region. Compared to the results presented in this 
paper, these applications require the calculation of further production vertices; 
however, for the treatment of the different scales and of the unintegrated gluon 
density all basic ingredients have been collected in this paper. 

Acknowledgements: A.S.V. thanks the Alexander-von-Humboldt Founda- 
tion for financial support. F.S. is supported by the Graduiertenkolleg "Zukiinftige 
Entwicklungen in der Tcilchcnphysik" . Helpful discussions with V. S. Fadin and 
L. N. Lipatov are gratefully acknowledged. 

A Alternative subtraction term 

In this Appendix we present an alternative subtraction term which does not 
make use of the simplifications A( 3 ) + Au\ — > 2Au\ and Ai$\ + Ar 6 ) — > 2A( 5 ) 
which we used in Eqs. (|63l l65f . These limits are valid in the kinematic regions 
leading to IR-divergences and hence they do provide the correct e poles. How- 
ever, they also alter the finite terms. Here we want to study also this finite 
part as accurately as possible and hence we do not take these limits but use the 
complete sum 

-4(i) + A ( 2) + A( 3 ) + A (4) + A {5) + A(fi) + A M rk (105) 
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as the gluonic subtraction term. 

The full gluonic matrix element written in Eq. (|54p contains spurious UV- 
divergences which are cancelled when combined with the MRK contribution. 
One fourth of the MRK contribution cancels the UV-divergence of ^4(4) while an- 
other fourth cancels that of A/gy The remaining half cancels the UV-divergence 
of two terms present in Eq. dSj) : 

A _ q'qh ( i-x i x i \ 

A(7) = " — {~rwi + r^kpj (106) 

which are IR-finite and hence so far not included in the subtraction term. 

By doubling Au\ and A(5) in the subtraction term constructed in Eq. (|67jl 
also their spurious UV-divergences are doubled and thus completely cancelled 
by the MRK contribution. But Eq. (|105[) so far only contains half of the spu- 
rious UV-divergences of the full matrix element in such a way that half of the 
MRK contribution is not compensated. Therefore a subtraction term based on 
Eq. (|105|) which is also free from spurious UV-divergences should also include 
A(j^ and A^ and reads 



= A {1) +A (2) +A (3) +A i4) +A i5) +A i6) + AmRK + A(7) + A( 8 ) 
= A {1) + A {2) + A m + (A {5) - A {4) ) + A {&) + + A {7) + A (8) . (108) 



If we now define 5(3,6,7,8) an d 5mrk as we did in Eq. (|88|) we get a new integrated 
subtraction term from the previous Eq. (|100p by replacing 
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l - (Shi + Siv) + S {3) + 5 (6) + + 5 (7) + S (8) . (110) 

The results for <S( 3 ) and <S( 6 ) can be easily obtained from Eqs. (C.43) and (C.40) 
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of Ref. 
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Due to the UV-singularity of Amrk we regularize the x integration by a cutoff 
8 to obtain 
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Making use of 2q a ki — q 2 = t + k 2 /x we can decompose Eq. (C.41) of Ref. [15] 
into one integration very similar to that of iSmrk an d another one which can be 
transformed to give 5(7) . 
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The two parts forming can be obtained from each other by the exchange 
ki <-> k% an d we only need to double the calculation of one: 
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When we add up these new contributions the spurious UV-divergences indeed 
cancel and we can safely take the <5 > limit. Furthermore, the new subtraction 
term has the same pole structure and only different finite parts when compared 
to that in Eq. ([67]) and its integrated form in Eq. (|100D . To complete the 
calculation we combine it with the corresponding unmodified quark part and 
obtain 
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